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Abstract 

We study the potential of the CNGS beam in constraining the parameter space of a 
model with one sterile neutrino separated from three active ones by an C(eV 2 ) mass- 
squared difference, Am^. We perform our analysis using the OPERA detector as 
a reference (our analysis can be upgraded including a detailed simulation of the 
ICARUS detector). We point out that the channel with the largest potential to 
constrain the sterile neutrino parameter space at the CNGS beam is i/n, — ► v T . The 
reason for that is twofold: first, the active-sterile mixing angle that governs this 
oscillation is the less constrained by present experiments; second, this is the signal 
for which both OPERA and ICARUS have been designed, and thus benefits from 
an extremely low background. In our analysis we also took into account i/„ — > v e 
oscillations. We find that the CNGS potential to look for sterile neutrinos is limited 
with nominal intensity of the beam, but it is significantly enhanced with a factor 2 
to 10 increase in the neutrino flux. Data from both channels allow us, in this case, 
to constrain further the four-neutrino model parameter space. Our results hold for 
any value of Am^ 0.1 eV 2 , i.e. when oscillations driven by this mass-squared 
difference are averaged. We have also checked that the bound on #13 that can be 
put at the CNGS is not affected by the possible existence of sterile neutrinos. 
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1 Introduction 



The results of solar PQEEEISin] , atmospheric [71IS] . reactor [9"|1U|I11|I12] and accelera- 
tor [T3lll4|ll5] neutrino experiments show that flavour mixing occurs not only in the 
hadronic sector, as it has been known for long, but in the leptonic sector as well. The 
full understanding of the leptonic mixing matrix constitutes, together with the dis- 
crimination of the Dirac/Majorana character of neutrinos and with the measurement 
of their absolute mass scale, the main goal of neutrino physics for the next decade. 

The experimental results point to two very distinct mass-squared differences, Am s 2 ol ~ 
7.9 x 1(T 5 eV 2 and |Am 2 tm | « 2.4 x 10~ 3 eV 2 . On the other hand, only two out of 
the four parameters of the three-family leptonic mixing matrix £/pmns [ToTll7|ll8|rT9"] are 
known: 6u ~ 34° and 823 ~ 43° [20]. The other two parameters, #13 and 5, are still 
unknown: for the mixing angle #13 direct searches at reactors [9~|10JI11] and three-family 
global analysis of the experimental data give the upper bound 6* 13 < 11.5°, whereas for 
the leptonic CP-violating phase 5 we have no information whatsoever (see, however, 
Ref. (23). 

The LSND data [2T|22II23] , on the other hand, would indicate a — > v e oscillation with 
a third neutrino mass-squared difference: AmL SND ~ 0.3 — 6 eV 2 , about two orders of 
magnitude larger than Am 2 tm . Given the strong hierarchy among the solar, atmospheric 
and LSND mass-squared splittings, Am 2 ol «C Am 2 tm «C Am 2 BL , it is not possible to 
explain all these data with just three massive neutrinos, as it has been shown with 
detailed calculations in Ref. [21] • A necessary condition to explain the whole ensemble 
of data in terms of neutrino oscillations is therefore the introduction of at least a 
fourth light neutrino state. This new light neutrino must be an electroweak singlet [18] 
in order to comply with the strong bounds on the Z° invisible decay width [251126] . For 
this reason, the LSND signal has often been considered as an evidence of the existence 
of a sterile neutrino. 

In recent years, global analyses of solar, atmospheric, short-baseline [27,28, 29, 30J exper- 
iments and LSND data have been performed to establish whether four-neutrino models 
can really reconcile the data and solve the puzzle [31ll32|33}34 .35.36. 37ll38j . The point 
is that providing a suitable mass-squared difference to each class of experiments is not 
enough: it is also necessary to show that the intrinsic structure of the neutrino mixing 
matrix is compatible with all the data. This turned out to be very hard to accomplish. 
In Ref. [32] it was shown that four-neutrino models were only marginally allowed, with 
best fit around AmLSND — 1 e V 2 and sin 2 26*lsnd — 1CT 3 . Generically speaking, the 
global analysis indicated that a single sterile neutrino state was not enough to reconcile 
LSND with the other experiments. For this reason, to improve the statistical compat- 
ibility between the LSND results and the rest of the oscillation data, models with 
two sterile neutrino states have been tested (see, for example, Ref. [10] and references 
therein). Although a slightly better global fit was achieved, a strong tension between 
the LSND data and the results from atmospheric and short-baseline experiments was 
still present. 



2 



So far, the LSND signal has not been confirmed by any other experiment |41j. It is 
therefore possible that the LSND anomaly arises from some some yet unknown problem 
in the data set itself. To close the issue, the MiniBooNE collaboration [42J at FermiLab 
has recently performed a search for — > v e appearance with a baseline of 540 m and a 
mean neutrino energy of about 700 MeV. The primary purpose of this experiment was 
to test the evidence for — > v e oscillation observed at LSND with a very similar L/E 
range. No evidence of the expected signal has been found, hence ruling out once and for 
all the four-neutrino interpretation of the LSND anomaly. However, MiniBooNE data 
are themselves not conclusive: although no evidence for — > v e oscillation has been 
reported in the spectrum region compatible with LSND results, an unexplained excess 
has been observed for lower energy neutrinos. Furthermore, within a five-neutrino model 
this excess can be easily explained, and even reconciled with LSND and all the other 
abearance experiments [33]. On the other hand, a post-MiniBooNE global analysis 
including also disappearance data show that five-neutrino models suffer from the same 
problems as four-neutrino schemes, and in particular they are now only marginally 
allowed - a situation very similar to that of four-neutrino models before MiniBooNE 
data. Adding a third sterile neutrincPI does not help [33] , and in general global analyses 
seem to indicate that sterile neutrinos alone are not enough to reconcile all the data. 
Models with sterile neutrinos and exotic physics have been therefore proposed (see, for 
example, Ref. [36]). 

In summary, the present experimental situation is still confused. It is therefore worth- 
wile to understand if, aside of MiniBooNE, new neutrino experiments currently running 
or under construction can investigate the existence of sterile neutrinos separated from 
the active ones by 0(eV 2 ) mass-squared differences. In this paper we explore in detail 
the capability of the CNGS beam to perform this search. For definiteness we focus on 
the simplest case with only one extra sterile neutrino. Note that this model is perfectly 
viable once the LSND result is dropped, as it contains as a limiting case the usual 
three-neutrino scenario. Furthermore, it is easily generalizable by adding new sterile 
neutrino states, and it can be used as a basis for models with extra "sterile" states 
strongly decoupled from active neutrinos (such as in extra-dimensions models with a 
right-handed neutrino in the bulk [37]). 

The CNGS beam j3H] has been built to test the (supposedly) dominant oscillation in at- 
mospheric neutrino data, — > v T . In order to make possible r production through CC 
interactions, the mean neutrino energy, (E v ) = 17 GeV, is much above the atmospheric 
oscillation peak for the CERN to Gran Sasso baseline, L = 732 Km. Two detectors are 
illuminated by the CNGS beam: OPERA (see Ref. [49J and refs. therein) will start data 
taking with the lead-emulsion target in 2007; ICARUS-T600 (see Ref. [SU] and refs. 
therein) will start operating in 2008. Both detectors have been especially designed to 
look for r's produced through —>■ v T oscillation and to minimize the corresponding 

1 A quite interesting scenario is, in our opinion, that in which three right-handed Majorana 
neutrinos are added to the three weakly-interacting ones. If the Majorana mass term M is 
O(eV), (3+3) light Majorana neutrinos are present at low-energy [44)45] . 
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backgrounds. The expected number of r events after signal selection in an experiment 
such as OPERA (after five years of data taking with nominal CNGS luminosity) is 
O(10) events with 0(1) background event. 

At the CNGS distance and energy, neutrino oscillations mediated by an 0(eV 2 ) mass 
difference will appear as a constant term in the oscillation probability. In four-neutrino 
models, fluctuations induced by this term over the atmospheric — > v T oscillation can 
be as large as 100% for specific points of the allowed parameter space. This is due to the 
fact that the leading angle for this oscillation is the less constrained one. The — > v T 
channel, therefore, is extremely promising as a "sterile neutrino" smoking gun, as it has 
been commented elsewhere (see, for example, Refs. [5Tf52] and refs. therein). To test 
the model we will also make use of the — > u e channel. Notice that the background to 
this signal coming from r — > e decay is modified in four-neutrino models with respect 
to standard three-family oscillations. In fact, since — > v T oscillations are depleted by 
active-sterile mixing with respect to standard ones, the r — > e background to — > v e 
oscillations gets depleted, too. A combined analysis of the two channels in four-neutrino 
models at the OPERA detector has been performed, taking into account properly all 
of the backgrounds. We stress, however, that the same analysis could be performed 
at ICARUS, as well. The previous considerations hold for any facility operating well 
beyond the kinematical threshold for r production. 

In the specific case of the CNGS beam, the limited flux implies a modest improvement 
in the parameter space exclusion, see Sec. El An increase in the exposure of such 
facilities, however, would permit to improve the present bounds on the parameters of 
four-neutrino models and, in particular, to constrain the leading angle in — > v T 
oscillations at the level of the other mixing parameters. 

The paper is organized as follows. In Sec. [2] we briefly review the main features of 
four-neutrino models and we introduce our parametrization of the mixing matrix. In 
Sec. [3] we compute the vacuum oscillation probabilities in the atmospheric regime and 
we review the present bounds on the active-sterile mixing angles. In Sec. H] we recall the 
most relevant parameters of CNGS. In Sec. [5] we study theoretically the expectations 
of the — > v T and — > v e channels at the CNGS. In Sec. |6] we present our results 
using these channels at the OPERA detector and the CNGS beam. Finally, in Sec. [7] 
we draw our conclusions. 



2 Four neutrino mass schemes 

In four-neutrino models, one extra sterile state is added to the three weakly interacting 
ones. The relation between the flavor and the mass eigenstates is then described by a 
4x4 unitary matrix U, which generalizes the usual 3x3 matrix ?7pmns [Tfill7|18fT5] . As 
stated in the introduction, in this work we only consider the case when the fourth mass 
eigenstate is separated by the other three by an 0(eV 2 ) mass-squared gap. There are 
six possible four-neutrino schemes, shown in Fig. (TJ that can accommodate the results 
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Fig. 1. The two classes of four-neutrino mass spectra, (3+1) and (2+2). 

from solar and atmospheric neutrino experiments and contain a third much larger Am 2 . 
They can be divided in two classes: (3+1) and (2+2). In the (3+1) schemes, there is 
a group of three close-by neutrino masses that is separated from the fourth one by 
the larger gap. In (2+2) schemes, there are two pairs of close masses separated by the 
large gap. While different schemes within the same class are presently indistinguishable, 
schemes belonging to different classes lead to very different phenomenological scenarios. 

A characteristic feature of (2+2) schemes is that the extra sterile state cannot be 
simultaneously decoupled from both solar and atmospheric oscillations. To understand 
why, let us define 



Vs 



E 

i £ sol 



and 



E 1^1 

j G atm 



(1) 



where the sums in i and j run over mass eigenstates involved in solar and atmospheric 
neutrino oscillations, respectively. Clearly, the quantities r) s and c s describe the fraction 
of sterile neutrino relevant for each class of experiment. Results from atmospheric and 
solar neutrino data imply that in both kind of experiments oscillation takes place 
mainly between active neutrinos. Specifically, from Fig. 46 of Ref. [20J we get rj s < 0.30 
and c s < 0.36 at the 3a level. However, in (2+2) schemes unitarity implies i] s + c s = 1, 
as can be easily understood by looking at Fig. [TJ These models are therefore ruled out 
at a very high confidence level |53j, and in the rest of this work we will not consider 
them anymore. 

On the other hand, (3+1) schemes are not affected by this problem. Although the 
experimental bounds on i] s and c s quoted above still hold, the condition i] s + c s = 1 no 
longer applies. For what concerns neutrino oscillations, (3+1) models are essentially 
unfalsifiable, since they reduce to the conventional three-neutrino scenario when the 
mixing between active and sterile states are small enough. 

The mixing matrix U can be conveniently parametrized in terms of six independent 
rotation angles dij and three (if neutrinos are Dirac fermions) or six (if neutrinos are 
Majorana fermions) phases 5i. In oscillation experiments, only the so-called "Dirac 
phases" can be measured, the effect of the "Majorana phases" being suppressed by 
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factors of m u /E u . The Majorana or Dirac nature of neutrinos can thus be tested only 
in AL = 2 transitions such as neutrino-less double /3-decay [51] or lepton number 
violating decays . In the following analysis, with no loss in generality, we will restrict 
ourselves to the case of 4 Dirac-type neutrinos only. 

A generic rotation in a four-dimensional space can be obtained by performing six dif- 
ferent rotations along the Euler axes. Since the ordering of the rotation matrices R^ 
(where ij refers to the plane in which the rotation takes place) is arbitrary, plenty of 
different parametrizations of the mixing matrix U are allowed. The large parameter 
space (6 angles and 3 phases, to be compared with the standard three-family mixing 
case of 3 angles and 1 phase) is however reduced to a subspace whenever some of 
the mass differences become negligible. If the eigenstates i and j are degenerate, ro- 
tations in the zj-plane become unphysical and the corresponding mixing angle should 
drop from oscillation probabilities. If the matrix is the rightmost one the angle 
8y automatically disappears, since the matrix commutes with the vacuum hamilto- 
nian. The parameter space gets therefore reduced to the physical angles and phases. 
If a different ordering of the rotation matrices is taken, no angle explicitly disappears 
from the oscillation formulas, but the physical parameter space is still smaller than the 
original one. In this parameter redefinition is needed to reduce the parameter 

space to the observable sector. In Refs. [S5.56J it was shown how the one-mass domi- 
nance (Agoi — > and A atm — > 0, where A = Am 2 L/AE |57j) and two-mass dominance 
(A so i — > 0) approximations can be implemented in a transparent way (in the sense that 
only the physical parameters appear in oscillation probabilities) using a parametriza- 
tion in which rotations are performed in the planes corresponding to smallest mass 
difference first: 



UsBL = Ru(0u) #24(024) -^34(6*34) #23(023, £3) #13(0133 £2) #12(012, #l) • (2) 



This parametrization was shown to be particularly useful when maximizing oscilla- 
tions driven by a 0(eV 2 ) mass difference. The analytical expressions for the oscillation 
probabilities in the (3+1) model in the one-mass dominance approximation in this 
parametrization have been presented in Ref. [5T] . 




In this paper, however, we are interested in a totally different regime: the "atmospheric 
regime", with oscillations driven by the atmospheric mass difference, Am^ tm L/£ ~ 
ii / 2. We will then make use of the following parametrization, adopted in Ref. [13] : 



U a tm = #34(034) #24(024) #23(023, £3) #14(014) #13(013, <5 2 ) #12(012, #l) • (3) 



It is convenient to put phases in i? 12 (so that it automatically drops in the two-mass 
dominance regime) and R13 (so that it reduces to the "standard" three-family Dirac 
phase when sterile neutrinos are decoupled). The third phase can be put anywhere; 
we will place it in R23. Note that in the one- mass dominance regime all the phases 
disappear from the oscillation probabilities. 
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Fig. 2. Allowed regions at 90%, 95%, 99% and 3a CL in the (0 13 , 14 ) plane (left) and in the 
(#24> #34) plane (right) from the results of present atmospheric, reactor and LBL neutrino 
experiments. The undisplayed parameters #23 and £3 are marginalized. 

3 Oscillation probabilities and allowed parameter space 

Let us first consider v e disappearance at L/E such that A so i can be safely neglected 
with respect to A atm and A SBL . We get for this probability in vacuum: 



sin 2 2#i4 sin^ 



AortX 



— c 14 sin 2 26*13 sin' 



A 



atm J 



(4) 



where c 



cos 6ij and 



sin 6ij. It is clear from Eq. (jlj) that reactor experiments 
such as Bugey and Chooz can put stringent bounds to 6*i 3 and 6 ) 14 , in this parametriza- 
tion. This is depicted in Fig. [2^1eft), where 90%, 95%, 99% and 3<r CL contours in 
the (6*i3 — 6 ) i 4 )-plane are shown for A so i — > and Am 2 tm = 2.4 x 10~ 3 eV 2 . The third 
mass difference, Am 2 BL , is free to vary above 0.1 eV 2 . Notice that the u e disappearance 
probability does not depend on #23, #24 and #34 . It can be clearly seen that the three- 
family Chooz bound on # 13 is slightly modulated by #14. Both angles, however, cannot 
be much larger than 10°. We will therefore expand in these two parameters to deduce 
the other relevant oscillation probabilities. 

At the CNGS beam atmospheric oscillations are large, solar oscillations can be ne- 
glected and C(eV 2 ) oscillations are extremely fast and can be averaged. It is useful to 
write down the oscillation probability (in vacuum) at typical atmospheric L/E, in the 
approximation A so i — > 0, A SBL — > 00. In this regime: 
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where + stands for neutrinos and — for antineutrinos, respectively. Up to second order 
in #13 and 6*14 we get for the disappearance oscillation probability: 



P. 



1 2c 14 s 2 4(l 
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S 23 C 24 



C 24l C 23 S 13; 
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2c2 4 s 23 (l - 2s2 3 )si 3 Si4S24 cos(<5 2 - 5 3 )} sin ; 
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atm J 



(6) 



A "negative" result in an atmospheric L/E disappearance experiment (such as, for 
example, K2K), in which oscillations can be very well fitted in terms of three-family 
oscillations, will put a stringent bound on the mixing angle #24 . The bound from such 
experiments on #24 can be seen in Fig. fright), where 90%, 95%, 99% and 3a CL 
contours in the (#24 — #34)-plane are shown for A so i — > and Am 2 tm = 2.4 x 10~ 3 eV 2 . 
The third mass difference, Am 2 BL , is free to vary above 0.1 eV 2 . The mixing angles not 
shown have been fixed to: #23 = 45°; 6*13 = #14 = (in this hypothesis, ? w does not 
depend on phases). Notice that the disappearance probability does not depend on 

#34- 

From the figure, we can see that #24 cannot be much larger than 10°, either. We will 
consider, therefore, the three mixing angles #13, #14 and #24 being of the same order and 
expand in powers of the three. At second order in #13, #14 and #24, we get: 



P — 1 — 2s 2 — 4s 2 



c 23 (l 2s 24 ) s 13 



. 2 ^atm-^ fr7 \ 
sin . (7) 



Since both v e and disappearance do not depend on #34, we should ask which mea- 
surements give the upper bound to this angle that can be observed in Fig. fright). 
This is indeed a result of indirect searches for — > v s conversion in atmospheric exper- 
iments, using the different interaction with matter of active and sterile neutrinos. This 
can be understood from the (vacuum) — > v s oscillation probability at atmospheric 
L/E for which, at second order in #13, #14 and #24, we get: 



5 D 



P - Or 2 c 2 _L_ /«ln 2 9fl (r A r 2 « 2 - r 2 e 2 
1 tis — z,c 34 15 24 ' \ alLL z ' L/ 23l 13°24' 5 34 "-"34*5 

+ 2c 34 sin 2^23^34 s 24 (l - 2s2 3 ) cos5 3 + 2s 23 Si3Si4COS(52 I sin 2 ^ a ^ mJ 
± c 34 sin 26 , 23 S24S34 sin S 3 sin A atm L . 



As it can be seen, the bound on #34 arises from a measurement of spectral distortion 
(i.e., from the "atmospheric" term proportional to sin 2 A atm L/2). On the other hand, 
bounds on #13, #14 and #24 are mainly drawn by a flux normalization measurement. 
As a consequence, the bound on #34 that we can draw by non-observation of — > v s 
oscillation in atmospheric experiments is less stringent than those we have shown before. 
For this reason, #34 can be somewhat larger than #13, #14 and #24, but still bounded to 
be below 40°. In the following, we will expand in powers of the four mixing angles 
#13, #14, #24 and #34, that will be considered to be comparably small. 
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Up to fourth-order in 6*13, 6*14, 624 and 6*34, the — > v e appearance probability in the 
atmospheric regime is: 

P?e = 4 {s2 3 Si 3 [l - s\ 3 - s 2 u - s|J + S23S13S14S24 cos(5 2 - <5 3 )} sin 2 Aat ™ L 
± 2S23S13S14S24 sin(5 2 - S 3 ) sin A atm L + 2s 2 4 S2 4 . 

Eventually, the — > u T appearance probability up to fourth-order in #13, #14, 6*24 and 
#34 in the atmospheric regime is: 



P^r — 2^24^34 + I sin 2 26> 2 3 [c\ 3 c\ A c 2 M — 

- 4 sin 26 l 23Si 3 si4 [s 23 s 3 4 cos 5 2 + c 23 s 2 4 cos(5 2 - S3)) 

A T (10) 

iO-oZ) 22 r 2 22i r 1 - 2 ^atm-^ v ' 

+ 2 sin 26> 2 3S24S34C 13 c 24 C34 |c 14 - 2c 13 s 23 J cosd 3 j sin — - — 
=F sin 26*23 S24S34C13C14C24C34 sin 5 3 sin A atm L . 

As it was shown in Refs. [SHE2], the — > v T appearance channel is a good place to 
look for sterile neutrinos. This can be understood as follows: consider the z/ M — > v T 
three-family oscillation probability in the atmospheric regime, up to fourth-order in 
#13 : 

A atir ,L 

T 

The genuine active-sterile neutrino mixing effects are: 



P% = P„ T (0 i4 = 0) ~ 4 sin 2 2023 sin 2 . (11) 



AP = P - p 3u 

L -^ 1 fir — 1 iir 1 MT 

= {-(4 + 4) sin 2 2^23 + 2 sin 2023(1 - 2s 2 3 ) S 24S34 cos<5 3 } sin 2 (12) 
=F sin 2023S24S34 sin 5 3 sin A atm L + . . . 

that is second-order in small angles 9 13 , 6*14, 6*24 and 6*34. We would get a similar result 
for Up disappearance, also. On the other hand, computing the corresponding quantity 
in the — > v e channel, we get: 

APfie = Pfie Pfj,e{@i4 = 0) 

= S23S13S14S24 cos(<5 2 - S 3 ) sin 2 a in (13) 
± 2S23S13S14S24 sin(5 2 - 5 3 ) sin A atm L + . . . 

that is third-order in the same parameters. 

Notice, eventually, that all oscillation probabilities start with an energy-independent 
term and are, therefore, non-vanishing for L = 0, a result of our assumption that 
A SRT , — >■ oo. 
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Fig. 3. CNGS neutrino fluxes (in arbitrary units) as a function of the neutrino energy. Both 
muon and electron neutrino fluxes are illustrated. 



4 The CNGS facility 



The CNGS is a conventional neutrino beam in which neutrinos are produced by the 
decay of secondary pions and kaons, obtained from collisions of 400 GeV protons from 
the CERN-SPS onto a graphite target. The resulting neutrinos are aimed to the under- 
ground Gran Sasso Laboratory (LNGS), located at 730 Km from CERN. This facility 
provided the first neutrinos in August 2006 [49]. Differently from other long baseline 
experiments, the neutrinos from CNGS can be exploited to search directly for — > v T 
oscillations, since they have a mean energy well beyond the kinematic threshold for r 
production. Moreover, the prompt v T contamination (mainly from D s decays) is negligi- 
ble. The expected v e contamination is also relatively small compared to the dominant 
component, thus allowing for the search of sub-dominant — > u e oscillations through 
an excess of v e CC events. 

The energy spectra of the CNGS neutrino beam are shown (in arbitrary units) in 
Fig. [3] [5H]. In the present paper we assume the nominal intensity for the CNGS, 
corresponding to 4.5 x 10 19 pot/year. 

OPERA has been designed to search for r appearance through identification of the 
v T CC interaction on an event-by-event basis. In particular, r's are tagged identifying 
explicitly their decay kink through high resolution nuclear emulsions interleaved with 
lead sheets. For this detector, we can take advantage of the detailed studies of the 
— > u T signal (see Ref. [59]) and of the — > v e signal (see Ref. |60j). 

The total non-oscillated CC event rates for a 1 Kton lead target with a neutrino flux 
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669.0 


13.7 


5.9 


0.3 



Table 1 

Nominal performance of the CNGS reference beam \58\j . The total non-oscillated CC event 
rates are calculated assuming 10 19 pot and 1 Kton lead target mass. 

normalized to 10 19 pot are shown in Tab. [T] and are evaluated according to 

N = j d -^^(E)dE, (14) 

in which <p Ua is the flux of the neutrino flavour u a and a Va the corresponding cross 
section on lead. 



5 Appearance channels at the CNGS 



5.1 — > v T oscillations 

Since the CNGS experiments have been designed to search for — > v T oscillation 
in the parameter region indicated by the atmospheric neutrino data, we can take full 
advantage of them in order to constrain (and, possibly, study) the four-family parameter 
space. 

The number of taus from — > v T oscillations is given by the convolution of the 
flux d(j) v {E)j dE with the v T charged-current cross-section on lead, a„ (E), weighted 
by the —>■ v T oscillation probability, P I1T (E), times the efficiency for the OPERA 
detector, e^ T : 

N, T = Aj P» T (E)a™(E) e, T dE . (15) 

A is a normalization factor which takes into account the target mass and the nor- 
malization of the Up flux in physical units. Specializing our analysis for the OPERA 
detector, we have considered an overall efficiency ~ 13%, [59]. This efficiency takes 
into account that OPERA is able to exploit several decay modes of the final state r, 
using both so-called short and long decays. 

The dominant sources of background for the —>■ v T signal are charm decays and 
hadronic reinteractions. Both of them only depend on the total neutrino flux and not 
on the oscillation probabilities. The OPERA experiment at the CNGS beam has been 
designed precisely to measure this channel, and thus the corresponding backgrounds 
are extremely low. 

In Tab. |2]we report the expected number of r events in the OPERA detector, according 
to Eq. (fl5l) . for different values of #13, #14, #24 and 6*34. Input points have been chosen 
according to the allowed regions in the parameter space shown in Sec. [3j The other 
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background 


(#13; #14! #24; #34) 


N T 


background 


(5°; 5°; 5°; 20°) 


8.9 


1.0 


(10°; 5°; 5°; 20°) 


8.5 


1.0 


(5°; 5°; 5°; 30°) 


6.9 


1.0 


(10°; 5°; 5°; 30°) 


6.5 


1.0 


(5°; 5°; 10°; 20°) 


8.3 


1.0 


(10°; 5°; 10°; 20°) 


7.9 


1.0 


(5°; 5°; 10°; 30°) 


10.5 


1.0 


(10°; 5°; 10°; 30°) 


10.3 


1.0 


3 families 


15.1 


1.0 


3 families 


14.4 


1.0 



Table 2 

Event rates and expected background for the — > v T channel in the OPERA detector, for 
different values of #14, #24 and #34 in the (3+1) scheme. The other unknown angle, #13 has 
been hxed to: #13 = 5°, 10°. Tiie CP-violating phases are: 81 = 82 = 0; £3 = 90°. As a 
reference, the expected value in the case of standard three-family oscillation (i.e., for 9^ = 0) 
is shown for maximal CP-violating phase 5. The rates are computed according to Eq. (|15p . 

parameters are: 12 = 34°; 8 23 = 45°; Am s 2 ol = 7.9 x lCT 5 eV 2 ; Am 2 tm = 2.4 x lCT 3 eV 2 
and Am 2 BL = 1 eV 2 (all mass differences are taken to be positive). Eventually, phases 
have been fixed to: 8\ = 82 — 0; £3 = 90°. The expected background is also shown. Rates 
refer to a flux normalized to 4.5 x 10 19 pot/year (the nominal intensity of the CNGS), 
an active lead target mass of 1.8 Kton and 5 years of data taking. For comparison, we 
also report the expected number of events in the usual 3-family scenario. 

As it can be seen, in most part of the parameter space we expect a significant depletion 
of the signal with respect to standard three-neutrino oscillations. However, the differ- 
ence between (3+1) model —>■ v T oscillations and standard ones is much bigger than 
the expected background. A good signal/noise separation can therefore be used to test 
the model. 

5.2 v^—* u e oscillations 

The number of electrons from the — > v e oscillation is given by the convolution of 
the flux d(j) v {E)/dE with the u e charged-current cross-section on lead, a^ e c (E), 
weighted by the —>■ u e oscillation probability, P^ e (E), times the efficiency for the 
OPERA detector, e^(E) [60]: 

N ^ = A j d -^r P ^( E X C ( E ) £ ^( E ) dE > ( 16 ) 

where A is defined as above. The overall signal efficiency is the convolution of the 
kinematic efficiency e^J 1 (that ranges from 60% to 80% for neutrino energies between 
5 to 20 GeV) and several (nearly factorizable) contributions. Among them, the most 
relevant are trigger efficiencies, effects due to fiducial volume cuts, vertex and brick 
finding efficiencies and the electron identification capability. They result in a global 
constant factor e{,": ct ~ 48%. 
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(013! 014! #24! #34) 


N e 




< C 


r — » e 




(5°; 5°; 5°; 20°) 


3.5 


19.4 


5.3 


2.8 


0.9 


(5°; 5°; 5°; 30°) 


3.5 


19.4 


5.3 


2.1 


0.9 


(5°; 5°; 10°; 20°) 


2.4 


19.4 


5.3 


2.3 


0.9 


(5°; 5°; 10°; 30°) 


2.4 


19.4 


5.3 


2.4 


0.9 


3 families 


3.7 


19.7 


5.3 


4.6 


0.9 


(10°; 5°; 5°; 20°) 


10.6 


19.4 


5.3 


2.7 


0.9 


(10°; 5°; 5°; 30°) 


10.4 


19.4 


5.3 


2.0 


0.9 


(10°; 5°; 10°; 20°) 


8.8 


19.4 


5.3 


2.2 


0.9 


(10°; 5°; 10°; 30°) 


8.6 


19.4 


5.3 


2.4 


0.9 


3 families 


15.1 


19.7 


5.3 


4.8 


0.9 



Table 3 

Event rates and expected background for the Un — > v e channel in the OPERA detector, for 
different values of #14, #24 and #34 in the (3+1) scheme. The other unknown angle, 9i 3 , has 
been fixed to: 6 13 = 5°, 10°. The CP-violating phases are: 5i = 5 2 = 0; 63 = 90°. As a 
reference, the expected value in the case of standard three-family oscillation(i.e., for On = 0) 
is shown for maximal CP-violating phase 5. The rates are computed according to Eq. (|16p . 
Backgrounds have been computed following Ref. [60]. 

The dominant sources of background to the — > v e signal are, in order of importance: 

(1) v e beam contamination; 

(2) fake electrons due to tc° decays from NC interactions; 

(3) electrons produced through r decay, where the r comes from z/ M — > v T oscillations; 

(4) CC Vfi events where the muon is lost and a track mimics an electron. 

Backgrounds (1), (2) and (4) depend very little on the oscillation parameters. On the 
other hand, the r — > e background depends strongly on the active-sterile mixing angles. 
As we have seen in Sec. 15.14 in the allowed region of the parameter space — > v T 
oscillations are significantly depleted with respect to the standard three-neutrino ones. 
As a consequence, this background gets depleted, too. 

In Tab.EJwe report the expected number of electrons in the OPERA detector, according 
to Eq. (1TB]) . for different values of #13, #14, 6*24 and 6*34. Input points have been chosen 
according to the allowed regions in the parameter space shown in Sec. [3j The other 
parameters are: 9 12 = 34°; 6 23 = 45°; Am s 2 ol = 7.9 x 1(T 5 eV 2 ; Am 2 tm = 2.4 x 1(T 3 eV 2 
and Am 2 BL = 1 eV 2 (all mass differences are taken to be positive). Eventually, phases 
have been fixed to: 8% = 5 2 = 0; 5 3 = 90°. Backgrounds have been computed accordingly 
to Ref. [60]. Rates refer to a flux normalized to 4.5 x 10 19 pot/year (the nominal intensity 
of the CNGS), an active lead target mass of 1.8 Kton and 5 years of data taking. For 
comparison, we also report the expected number of events in the usual 3-family scenario. 
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As it can be seen from Tab. [3j the difference between the (3+1) model and the stan- 
dard three-neutrino oscillations are smaller in this channel than in the — > v T one. 
Moreover, they linearly depends on 9 13 , as it is clear from Eq. ({TBI . For # 13 = 5°, this 
channel will be of no help to test the allowed parameter space of the (3+1) model. On 
the other hand, for #13 saturating the Chooz-Bugey bound, both — > u T and — > u e 
might cooperate. However, notice that backgrounds to this signal are much larger than 
the difference between (3+1) model and standard three-neutrino oscillations for any 
value of 6> 13 . 



6 Sensitivity to (3 + 1) sterile neutrinos at OPERA 

In this section we study the sensitivity to #13 and to the active-sterile mixing angles 
#14, #24 and 6*34 at the CNGS beam, using both the — > u T and — > u e appearance 
channels at the OPERA detector. In the rest of this section, the known three-family 
subspace angles have been fixed to: #12 = 34°; #23 — 45°. The mass differences have 
been fixed to: Am s 2 ol = 7.9 x 10" 5 eV 2 and Am 2 tm = 2.4 x 10~ 3 eV 2 . The CP- violating 
phases 5\ and 62 have been kept fixed to 5i = 62 = 0. On the contrary, the CP- violating 
phase £3 is fixed to two values: £3 = or 90°. Notice that this phase is still present in 
the oscillation probabilities even when # 12 and #13 vanish, see Eq. ffTUj) . At atmospheric 
L/E, oscillations driven by an 0(eV 2 ) mass difference are averaged. We have checked 
that our results apply for any value of Am 2 BL > 0.1 eV 2 . 

In Fig. H] we show the sensitivity limit at 99% CL in the (#13, #14) plane (left) and 
in the (#24, #34) plane (right) from a null result of the OPERA experiment, assuming 
1, 2, 3, 5 and 10 times the nominal intensity of 4.5 x 10 19 pot/year. The coloured 
regions show the present bounds at 90% and 99% CL. We assume #23 = 45° and 
63 = 0° (top) or 5 3 = 90° (bottom). The sensitivity is defined as the region for which a 
(poissonian) 2 d.o.f.'s x 2 is compatible with a "null result" at the 99% CL. We refer to 
"null result" when #13 and the three active-sterile mixing angles, #14, #24 and #34 vanish 
simultaneously. Both — > v T and — > v e oscillations have been considered, with the 
corresponding backgrounds treated properly as in Sec. [5j An overall systematic error 
of 10% has been taken into account. 

In the left panels of Fig. H]we can see that OPERA can improve only a little the bound 
on 613 after 5 years of data taking working at nominal CNGS beam intensity, both 
for S3 = (top panel) or £3 = 90° (bottom panel). Increasing the nominal intensity, 
however, a significant improvement on the bound is achieved for any value of #14. Notice 
that the limit on #14 is almost unaffected by the OPERA data. This is because for the 
v ii v t an d Vp, — > v e oscillation probabilities at atmospheric L/E, the ^-dependence 
always arises at third-order in the small parameters #13, 6*i 4 , # 2 4 and #34 (see Eqs. 
and (JTOl) for the explicit expression in the adopted parametrization, Eq. ([3])). On the 
contrary, the #13-, #24- and ^-dependences in the same oscillation probabilities are 
quadratic in the small parameters. In case of vanishing active-sterile mixing angles, 
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Fig. 4. Sensitivity limit at 99% CL in the (#13, #14) plane (left) and in the (824, 634) plane 
(right) from a null result of the OPERA experiment, assuming 1, 2, 3, 5 and 10 times the 
nominal intensity of 4.5 x 10 19 pot/year. The coloured regions show the present bounds at 
90% and 99% CL. We assume 6> 23 = 45° and 5 3 = 0° (top) or 5 3 = 90° (bottom). 

8 H = 0, see Ref. [60]. 

In the right panels of Fig. H] the sensitivity of OPERA to 824 and #34 is shown. First 
of all, notice that the sensitivity is strongly affected by the intensity of the beam. 
No improvement on the existing bounds on these two parameters is achieved after 5 
years of data taking at nominal CNGS beam intensity, for any of the considered value 
of 5 3 . Already with a doubled flux intensity, some sensitivity to 824, 834 is achievable. 
The sensitivity enhancement strongly depends on the value of the CP- violating phase 
63, however. For 5 3 = 0, OPERA can exclude a small part of the 99% CL allowed 
region, only. On the other hand, for <5 3 = 90° twice the nominal CNGS flux suffices 
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Fig. 5. Sensitivity limit at 99% CL in the (#13, #14) plane (left) and in the (#24, #34) plane 
(right) from the combined analysis of present data and a null result of the OPERA experiment, 
assuming 1, 2, 3, 5 and 10 times the nominal intensity of 4.5 x 10 19 pot/year. The coloured 
regions show the present bounds at 90% and 99% CL. We assume 823 = 45° and 63 = 0° 
(top) or 5 3 = 90° (bottom). 

to put a bound on 834 < 25° for #24 > 4° at 99% CL. For maximal CP-violating 
S3, increasing further the CNGS flux can significantly constrain the (624, O34) allowed 
parameter space. Notice, eventually, the strong correlation between #24 and 634 in the 
right panels of Fig. HI This is an indication that the dominant channel that constrains 
these angles is — > v T . As it can be seen in Eq. f JTUj) . the two angles always appear in 
combination, with an approximate exchange symmetry #24 *-* #34- 

The allowed regions at 99% CL in the (6*13, #14) plane (left) and in the (624, O34) plane 
(right) from the combined analysis of present data and a null result of the OPERA 
experiment after 5 years of data taking (assuming 1, 2, 3, 5 and 10 times the nominal 
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CNGS intensity of 4.5 x 10 19 pot/year) are eventually shown in Fig. [51 The coloured 
regions refer to the present bounds at 90% and 99% CL, for 9 2 3 = 45° and 5 3 = 0° (top) 
or £3 = 90° (bottom). As it can be seen, the sensitivity of OPERA strongly benefits 
from the complementary information on the neutrino parameters provided by other 
experiments. In this case, even with the nominal beam intensity the extension of the 
allowed regions is reduced by a moderate but non-negligible amount. 



7 Conclusions 



The results of atmospheric, solar, accelerator and reactor neutrino experiments show 
that flavour mixing occurs not only in the quark sector, as it has been known for 
long, but also in the leptonic sector. Experimental data well fit into a three-family 
scenario. The existence of new "sterile" neutrino states with masses in the eV range is 
not excluded, however, provided that their couplings with active neutrinos are small 
enough. 

In this paper, we have tried to test the potential of the OPERA experiment at the 
CNGS beam to improve the present bounds on the parameters of the so-called four- 
neutrino models. The model, in which only one sterile neutrino is added to the three 
active ones responsible for solar and atmospheric oscillations, is the minimal extension 
of the standard three-family oscillation scenario. 

We have determined the presently allowed regions for all active-sterile mixing angles 
and studied the OPERA capability to constrain them further using both the — > v e 
and — > v T channels. We have performed our analysis using the OPERA detector as 
a reference. It can be extended including a detailed simulation of the ICARUS detector 
at the CNGS beam. 

Our conclusions are the following: if the OPERA detector is exposed to the nominal 
CNGS beam intensity, a null result can improve a bit the present bound on #13, but 
not those on the active-sterile mixing angles, #14 , #24 and #34. If the beam intensity is 
increased by a factor 2 or beyond, not only the sensitivity to #13 increases accordingly, 
but a significant sensitivity to 6*24 and 6*34 is achievable. The (#24, #34) sensitivity strongly 
depends on the value of the CP- violating phase #3, however, with stronger sensitivity for 
values of £3 approaching tt/2. Only a marginal improvement is achievable on the bound 
on 614, that should be constrained by high-intensitiy v e disappearance experiments. 

Notice that our results hold for any value of Am| BL > 0.1 eV 2 , i.e. in the region of L/E 
for which oscillations driven by this mass difference are effectively averaged. 
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